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1.  Introduction 

Recent  theoretical  and  experimental  studies  of  mul¬ 
tiple  scattering  and  multiple-field-of-view  (MFOV)  li¬ 
dar  detection  have  made  possible  the  retrieval  of 
cloud  droplet  size.1-6  Basically  and  as  originally  in¬ 
dicated  by  Eloranta6  and  Platt,7  the  enhanced  lidar 
signal  caused  by  multiple  scattering  contains  infor¬ 
mation  about  droplet  size.  Multiple-scattering 
events  are  always  present  in  clouds  and  are  usually 
measured  by  detection  of  the  cross-polarized  compo¬ 
nent  of  the  received  signals8-9  or  by  measurement  of 
the  lidar  returns  at  different  fields  of  view.10  In  our 
previous  theoretical  model1  we  used  the  total  (sum  of 
both  polarizations)  MFOV  lidar  signal,  which  limits 
the  analysis  to  cloud  depths  large  enough  for 
multiple-scattering  contributions  to  be  significantly 
above  the  lidar  pulse-to-pulse  signal  fluctuation. 
This  means  that  droplet-size  retrieval  was  not  possi¬ 
ble  at  small  depths  of  penetration  into  the  cloud.  In 
this  paper  we  show,  both  theoretically  and  experi¬ 
mentally,  that  the  use  of  secondary  polarization 
makes  possible  the  retrieval  of  the  droplet-size  den¬ 
sity  function  of  clouds  of  small  optical  depth. 

In  the  absence  of  multiple  scattering,  the  lidar  sig- 
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nal  consists  of  backscattered  light  at  exactly  180°, 
and  a  linearly  polarized  light  source  produces  no 
secondary-polarization  lidar  signal  because,  for 
spherical  droplets,  the  depolarization  ratio  is  zero  at 
180°.  To  illustrate  this,  Figs.  1  and  2  show  the 
phase  function  and  the  depolarization  ratio  of  lin¬ 
early  polarized  light  obtained  from  Mie  calculations 
for  a  type  C2  water-droplet  cloud11  and  a  submi¬ 
crometer  oil-droplet  suspension.  The  high  peak 
values  of  the  phase  function  near  0°  for  the  water- 
droplet  cloud  are  attributed  to  diffraction  and,  there¬ 
fore,  contain  particle-size  information.  Diffraction 
does  not  cause  depolarization,  as  is  clearly  shown  in 
Fig.  1,  but  the  multiple  internal  reflections  in  the 
water  droplets  cause  Strong  depolarization  at  large 
scattering  angles.  The  phase  function  and  the  de¬ 
polarization  ratio  of  the  oil-droplet  cloud  are  typical  of 
Rayleigh  scattering,  with  a  signature  close  to  a  dipole 
(phase  function  quasi-flat  and  a  zero-depolarization 
ratio  except  at  90°).  As  we  stated  above,  both  depo¬ 
larization  ratios  of  Figs.  1  and  2  rapidly  drop  to  zero 
at  exactly  180°. 

2.  Model 

Basically,  MFOV  lidar  measurements  are  not  local 
measurements.  They  depend  on  the  forward- 
scattered  radiation  by  all  aerosol  particles  preceding 
the  point  where  the  backscattering  occurs.  The 
droplet-size  distribution  function,  particularly  the 
number  density,  at  the  backscattering  point  has  an 
effect  on  the  amplitude  of  the  backscatter  signal  but 
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Fig.  1.  Phase  function  and  depolarization  of  a  C2-type  cloud. 


has  little  effect  on  its  angular  spreading,  as  the  back- 
scattering  phase  function  is  rather  flat  in  comparison 
with  the  forward  peak  from  140°  to  180°.  Figure  3 
illustrates  the  scattering  processes  that  lead  to  sec¬ 
ondary  polarization  in  the  presence  of  multiple  scat¬ 
tering  and  defines  the  geometry  of  the  problem.  To 
keep  the  mathematics  tractable  we  do  not  consider 
scattering  orders  higher  than  2  in  this  model,  which 
has  the  effect  of  limiting  the  application  of  the  model 
to  small  optical  depth.  However,  the  effect  of  higher 
scattering  orders  is  discussed  in  Section  5  below. 
The  cloud  is  located  at  distance  za  from  the  lidar,  and 
a  first  forward  scattering,  along  (3,  occurs  at  a  dis¬ 
tance  z  followed  by  a  backscattering,  along  ir  -  |3  +  0, 
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Fig.  2.  Phase  function  and  depolarization  of  a  fog-oil  cloud. 


in  the  plane  zc.  The  field  of  view  (FOV)  0  and  the 
scattering  angle  (3  are  related  as  follows: 


tan  (3  = 


- 


tan  0. 


(1) 


\ZC  -  Zj 

Given  a  cloud  with  extinction  coefficient  a(z),  we 
obtain  the  measured  received  lidar  power  S  in  the 
secondary  polarization  between  FOV’s  0J+1  and  Qj  by 
summing  all  the  scattering  contributions  between 
distances  za  and  zc  over  particle-size  density  distri¬ 
bution  q0(x),  where  x  is  the  particles’  diameter;  i.e., 


S(zc,  QJ+1  -  0,0  =  P0  exp[-2a(zc 
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where  I(z,  x,  (3,  n)  is  the  forward-scattering  phase 
function,  8(z,  x,  <(3>,  n)  is  the  depolarization  ratio  at 
average  backscattering  angle  tt  —  <(3)  +  0,  ((3)  = 
[(37+1(z)  +  |3y(z)]/2,  n  is  the  refractive  index,  P0  is  the 
laser  power,  ct  is  the  product  of  the  speed  of  light  and 
the  pulse  width,  p  is  the  backscattering  coefficient  at 
zc,  and  a(z)  is  the  scattering  coefficient.  The  num¬ 
ber  density  distribution  q0(x)  can  be  converted  to  the 
volume  density  distribution12  q3(x)  by  means  of  the 
equation 


q0(x)  = 


x  3q3(x ) 


x  3q3(x)dx 


(3) 


and  the  integral  over  x  in  Eq.  (2)  is  divided  into  M 
particle-size  intervals  with  q3(x)  set  equal  to  a  con¬ 
stant  in  each  interval.  We  thus  rewrite  Eq.  (2)  as 
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Fig.  3.  Scattering  processes  that  lead  to  secondary  polarization. 


Equation  (4)  is  rewritten  in  matrix  form  as 

S  =  Aqs,  (5) 

where  the  matrix  coefficients  are  given  by 


’^i+l 

1* 

Zq 
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hj 

X  ot(z)8(z,  x,  ((3),  7i)dzdx,  (6) 

where  C3  is  a  constant  with  respect  to  z,  x,  and  (3. 
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Prom  Eq.  (5)  and  (6)  the  particle-size  density  dis¬ 
tribution  (PSD)  qg  is  obtained  through  matrix  inver¬ 
sion.  However,  as  in  many  inverse  problems  in 
physics,  the  direct  inversion  of  Eq.  (5)  does  not  lead  to 
satisfactory  results.  It  is  necessary  to  use  a  con¬ 
strained  linear  inversion  technique.  Equation  (6) 
contains  the  unknown  range-dependent  scattering 
coefficient  a(z).  However,  because  the  system  is  not 
calibrated  and  the  constant  C3  is  unspecified,  the 
PSD  can  be  determined  only  in  relative  units. 
Therefore  we  need  only  the  relative  strength  of  a(z). 
Because  the  method  is  applicable  at  small  penetra¬ 
tion  depths,  we  can  take  a(z)  proportional  to  the  mea¬ 
sured  total  lidar  signal. 

Before  considering  calculations,  we  first  need  to 
define  the  depolarization  8(z,  x,  ((3),  n),  and  the 
forward-scattering  phase  function  I(z,  x,  fi,  n).  The 
depolarization  8(2,  x,  ((3),  n)  is  obtained13  by  Mie  the¬ 
ory  and  is  given  by 

8(z,  x,  <(3),  n) 

P2(0»cos2(Q))  -  2P3(Q»Cos«3»  +  ^(((3)) 
3P2«p»cos2«(3»  +  2P3«|3»cos«|3})  +  SP^O))  ’ 

(7) 

where  P1,  P2,  and  P3  are  the  scattered  Stokes  param¬ 
eters. 

For  droplets  that  are  significantly  larger  than  the 
wavelength,  the  forward-scattering  phase  function 
I(z,  x,  (3,  n)  can  be  obtained  by  use  of  Mie  theory  or  the 
Fraunhofer  diffraction  theory14  in  conjunction  with 
an  analytical  expression  for  the  geometrical  optics 
component.16-16  Note  that  the  agreement  between 
the  two  theories  is  quite  good  for  droplets  larger  than 
the  wavelength  and  for  scattering  angles  smaller 
than  45°.  In  Eq.  (2)  the  integral  over  (3  corresponds 
to  a  difference  of  encircled  energy.  This  difference  is 
maximum  in  the  region  where  the  encircled  energy 
function  increases  the  most  rapidly,14  and  the 
particle-size  ranges  in  Eq.  (2)  are  chosen  in  this  re¬ 
gion  to  yield  the  maximum  energy  difference  for  the 
diagonal  matrix  elements,1  which  helps  in  achieve¬ 
ment  of  a  robustness  that  naturally  helps  in  the  in¬ 
version  of  matrix  A.  The  Fraunhofer  diffraction 
theory  leads  to  an  easy-to-use  analytical  definition  of 
the  encircled  energy,  and  it  is  useful  in  defining  the 
particle-size  range.  However,  Mie  theory  is  used 
here  to  extend  the  applicability  of  the  inversion  tech¬ 
nique  to  droplets  smaller  than  the  wavelength. 

Figure  4  shows  the  calculated  [Eq.  (4)]  range- 
corrected  s-polarization  lidar  return  within  FOV’s  0j 
and  0I+1  as  a  function  0,+1  for  a  uniform  cloud  begin¬ 
ning  at  a  distance  of  95  m  from  the  lidar  and  with  a 
penetration  depth  of  6  m.  The  assumed  water- 
droplet  size  density  distributions  are  three  log¬ 
normal  distributions,  with  20,  10,  and  0.8  |xm  as 
geometric  mean  diameters  and  with  geometric  stan¬ 
dard  deviations  lncr  =  0.2,  0.2,  0.58,  respectively. 
The  larger  the  droplet,  the  narrower  the  diffraction 
peak  and  therefore  the  maximum  the  energy  col- 


Fig.  4.  Range-corrected  s-polarization  lidar  return  within  FOV’s 
6;  and  0i+1  as  a  function  of  FOV  0i+1.  The  curves  are  for  log¬ 
normal  droplet-size  distributions  of  several  mean  diameters  and 
logarithmic  standard  deviations. 


lected  at  small  FOV’s.  The  FOV  signature  is  dis¬ 
tinct  and  characterizes  the  droplet-size  density 
distribution.  Figure  5  shows  the  corresponding  sig¬ 
nals  integrated  over  FOV  0,-+1  as  a  function  of  0i+1. 
These  curves  give  the  secondary-polarization  energy 
contained  in  FOV  0;+1. 

3.  Inversion 

Because  the  direct  inversion  of  Eq.  (5J  does  not  lead  to 
satisfactory  results  it  is  necessary  to  use  a  con¬ 
strained  linear  inversion  technique,17-18  which  gives  a 
solution  of  the  form 

q  =  (ATA  +  yHfWP,  (8) 

where  7  is  a  Lagrangian  multiplier  and  H  (the  second 
difference)  is  a  matrix  that  we  have  used  to  establish 
a  relation  between  successive  values  of  qt . 


Fig.  5.  Range-corrected  s-polarization  lidar  return  within  FOV 
0i+1  as  a  function  of  FOV  6 £ 1 . 
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Fig.  6.  Ideal  lidar  signals  S'(0;+1)  -  jS(0;)  and  signals  with  20% 
white  noise  as  function  of  FOV  0j+1  for  a  cloud  made  from  log- 
normally  distributed  particles  with  median  diameters  of  10  and  0.8 
pm.  The  distance  to  the  lidar  is  95  m,  and  the  penetration  depth 
is  6  m. 


One  determines  the  Lagrangian  multiplier  of  Eq. 
(8)  by  performing  the  inversion  on  synthetic  MFOV 
lidar  signals  calculated  for  a  known  PSD.  Typically, 
the  range  of  y  is  within  the  range  of  the  product  of  the 
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Fig.  7.  Effect  of  y  on  the  recovery  of  the  initial  density  distribu¬ 
tion:  (a)  xm  —  10  pm,  (b)  xrn  =  0.8  pm.  sig,  Logarithmic  stan¬ 
dard  deviation. 


matrix  coefficients.  Figure  6  shows  the  lidar  signals 
S(0£+i)  —  S(0j)  as  functions  of  FOV0£+1  for  two  log¬ 
normal  distributions,  i.e., 


qs(x)  =  ■  ■  r~-  exp[-  (In  x  -  In  xm)2/a2],  (9) 

xcryzir 

with  xm  =  0.8,  10  |xm  and  ct  =  0.58,  0.2,  respectively. 
The  water-droplet  cloud  is  located  at  a  distance  za  of 
95  m  from  the  lidar.  The  target  position  is  at  101  m 
for  a  penetration  depth  inside  the  cloud  of  6  m.  The 
refractive  index  of  water  is  equal  to  (1.33,  0).  In  the 
same  figure,  the  same  S(0£+1)  -  S(Qt)  are  replotted 
with  20%  rms  white  noise. 

Figures  7(a)  and  7(b)  show  the  particle-size  volume 
density  distributions  recovered  for  various  y  in  the 
absence  of  noise  for  xm  =  10  |xm  and  xm  =  0.8  pm, 
respectively.  Figure  8  illustrates  the  effect  of  added 
noise.  A  y  of  10 - 3  appears  to  be  an  acceptable  value. 
The  recovered  distributions  obtained  for  the  0.8-pm 
particles  are  not  so  good  as  those  obtained  for  the 
larger,  10-pm  particles  in  the  absence  of  noise. 
Small  particles  are  characterized  by  large  scattering 
angles,  and  the  inaccuracy  of  the  inversion  is  likely 
caused  by  undersampling  of  large  scattering  angles. 
For  example,  less  than  10%  of  the  cloud’s  probed 
length  contributes  to  measurements  at  scattering  an- 
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Fig.  8.  Influence  of  noise  on  choice  of  y:  (a)  xm  -  10  pm,  (b)  xm  = 
0.8  pm.  sig,  logarithmic  standard  deviation. 
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gles  greater  than  45°  for  our  measurements  made  at 
a  6-m  penetration  depth. 

4.  Experimental  Validation 

MFOV  lidar  measurements  were  made  under  a  con¬ 
trolled  environment  with  a  22-m-long  aerosol  cham¬ 
ber.  Two  dissemination  systems,  distributed 
pneumatic  water  nozzles  and  a  MDG  Super  Max 
5000  fog-oil  generator,  were  used  to  generate  droplets 
of  fairly  different  sizes.  In  the  first  case  the  gener¬ 
ated  water  droplets  were  10-80  |xm  in  diameter,  in 
the  second  case  the  oil  droplets  had  submicrometer 
diameters.  The  aerosol  chamber  was  situated  95  m 
from  the  lidar.  W e  used  longitudinal  and  transverse 
transmissometers  to  measure  the  optical  depth,  and 
in  situ  particle  sizing  of  the  water  droplets  was  per¬ 
formed  in  the  middle  of  the  chamber  with  a  Malvern 
particle  sizer.  Figure  9  shows  a  typical  water- 
droplet  volume  density  distribution  such  as  mea¬ 
sured  with  the  Malvern  particle  sizer.  In  the  same 
figure  we  also  show  a  fog-oil-droplet  volume  density 
distribution  obtained  from  the  literature19  (a  log¬ 
normal  distribution  with  a  geometric  mean  volume 
diameter  of  0.8  |xm  and  a  geometric  standard  devia¬ 
tion  of  0.58).  Fog-oil  droplets  are  too  small  to  be 
measured  with  the  Malvern  instrument. 

The  MFOV  lidar  measurements  were  made  se¬ 
quentially  with  a  100-Hz  repetition-rate  Nd:YAG  la¬ 
ser  (Kigre  Workhorse  Model  1732)  synchronized  with 
a  rotating  aluminized  glass  disk  etched  at  the  periph¬ 
ery  with  32  irises  or  rings  of  different  apertures,  as 
illustrated  in  Fig.  10.  The  ring-mask  measurements 
have  the  advantage  of  being  less  sensitive  to  shot-to- 
shot  cloud  fluctuations.  However,  the  small  area  of 
the  inner  rings  make  the  individual  measurements 
fall  below  the  detection  threshold  quickly  with  in¬ 
creasing  optical  depth  and  thus  lead  to  invalid  cumu¬ 
lative  measurements.  Table  1  lists  the  individual 
rings’  inner  and  outer  diameters.  The  diameters  of 
the  iris  masks  are  identical  to  the  outer  diameters  of 
the  ring  masks.  Both  disks  are  mounted  in  the  im- 


Fig.  9.  Volume  density  distributions  of  a  generated  water-droplet 
cloud  measured  with  a  Malvern  Particle  Sizer  and  of  typical  pub¬ 
lished  values  for  fog-oil  clouds.17 


(b) 


Pig.  10.  Aluminized  glass  plate  with  32  etched  irises,  (a)  The 
order  of  the  irises  was  set  to  cover  the  maximum  FOV  difference 
within  1/4  of  a  turn,  (b)  The  order  of  the  rings  is  the  same  as  for 
the  irises.  _ _ 


age  plane  of  the  primary  optics,  as  illustrated  in  Fig. 
11.  The  primary  optics  consists  of  a  200-mm- 
diameter  off-axis  parabolic  mirror  with  a  focal  length 
of  760  mm.  The  position  of  the  image  plane  is  a 
function  of  the  focal  length  and  the  object  position, 
and  it  is  necessary  to  adjust  the  image  plane  position 
in  accordance  with  the  cloud  distance.  Secondary 
optics  and  a  polarizer  cube  are  used  to  split  and  re¬ 
image  the  two  polarization  components  on  3-mm- 
diameter  avalanche  photodiodes  (APD’s).  We 
calibrated  the  radial  responses  of  the  APD  detectors 
by  scanning  the  laser  spot  on  a  canvas  target.  The 
glass  disk  turns  at  a  stabilized  speed  of  3.125  rota¬ 
tions/  s,  and  an  electronic  delay  was  set  between  the 
synchronization  holes  and  the  laser  trigger.  In  this 
configuration  the  laser  is  slaved  to  the  FOV  control¬ 
ler.  FOVs  ranging  from  0.1  to  12  mrad  were  thus 
changed  at  a  repetition  rate  of  100  Hz.  Finally,  the 
characteristics  of  the  laser  beam  were  as  follows: 
2.5-cm  diameter  and  0.3-mrad  divergence  (50%  total 
energy),  linear  polarization  purity  of  1/500  ensured 
with  a  high-power  polarizer  cube,  an  outgoing  energy 
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Table  1.  Inner  and  Outer  Ring  Diameters  Etched  on  the  Rotating 
Glass  Disk 


Ring 

Number 

Inner  Ring 
Diameter  (pm) 

Outer  Ring 
Diameter  (pm) 

1 

.  _  0 

76 

2 

76 

89 

3 

89 

104 

4 

104 

121 

5 

121 

142 

6 

142 

165 

7 

165 

193 

8 

193 

226 

9 

226 

264 

10 

264 

308 

11 

.308 

360 

12 

360 

420 

13 

420 

491 

14 

491 

574 

15 

574 

670 

16 

670 

783 

17 

783 

915 

18 

915 

1069 

19 

1069 

1249 

20 

1249 

1458 

21 

1458 

1704 

22 

1704 

1990 

23 

1990 

2326 

24 

2326 

2716 

25 

2716 

3174 

26 

3174 

3707 

27 

3707 

4331 

28 

4331 

5059 

29 

-5059 

5911 

30 

5911 

6905 

31 

6905 

8064 

32 

8064 

9424 

per  pulse  in  the  atmosphere  of  60  mJ,  and  a  pulse 
width  of  12  ns. 

The  measurements  were  usually  performed  early 
in  the  morning  so  the  relative  humidity  was  high  and 
the  wind  mild.  Table  2  lists  four  trials  that  we  used 


Table  2.  Trial  List  Parameters 


Trial  Number 
and  Medium 

Penetration 
Depth  (m) 

Optical 

Depth 

FOV 

Type 

pvri0318,  water 

6 

0.4 

Ring 

pvh0305,  water 

6 

0.2 

Iris 

pvh0305,  water 

10 

0.4 

Iris 

fori0308,  oil 

6 

0.4 

Ring 

foh0308,  oil 

6 

0.4 

Iris 

to  demonstrate  the  PSD  lidar  polarization  technique. 
Figure  12  shows  the  water-droplet  trial  number 
pvri0318  principal  p-  and  secondary  s-polarization 
lidar  returns  as  functions  of  distance  for  two  total 
FOV’s.  For  this  representation  the  individual  ring 
signals  are  summed  to  0.75  and  12  mrad,  respec¬ 
tively.  In  Fig.  13  we  display  the  corresponding  sig¬ 
nals  for  oil-droplet  trial  number  foh0301.  Clearly  at 
the  beginning  of  each  cloud  the  p-polarization  mea¬ 
surement  values  are  2-3  orders  of  magnitude  higher 
than  the  s-polarization  values.  Without  polariza¬ 
tion  separation,  the  multiple-scattering  signal  would 
be  lost  in  the  main  polarization  signal  fluctuations. 
For  the  large  water  droplets,  the  differences  between 
the  two  FOV  signals  increase  with  distance  (optical 
depth),  whereas  for  the  small  oil  droplets  the  differ¬ 
ences  remain  quasi-constant,  as  is  particularly  evi¬ 
dent  for  s  polarization. 

Figures  14  and  15  show  the  measured  secondary 
s-polarization  ring  and  cumulative  signals,  respec¬ 
tively,  as  functions  of  FOV  e,+i.  at  a  distance  of  101  m 
(i.e.,  6  m  into  the  cloud)  for  the  water-  and  fog  oil 
droplet  clouds.  Figure  16  shows  the  cumulative 
measurements  obtained  for  trials  pvh0305  and 
foh0301  at  distances  of  101  and  105  m  (i.e.,  6  and 
10  m  into  the  cloud).  The  profiles  of  Figs.  14-16  are 
multiple-scattering  measurements  and  show  trends 
comparable  to  those  of  the  model-based  calculated 
signals  plotted  in  Figs.  4  and  5.  To  better  illustrate 
the  effect  of  multiple  scattering,  we  added  to  Fig.  15 


Laser  line 


Fig.  11.  Setup  for  the  detection  module  used  for  sequential  MFOV  lidar  measurements.  FOV  control  is  basically  the  ring  or  iris  rotating 
dish  assembly  mounted  on  a  precision  translator. 
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Pig.  12.  Range-corrected  lidar  return  p-  and  secondary 
s-polarization  lidar  return  as  functions  of  distance  for  0.75-  and 
12-mrad  total  FOVs  for  water  droplet  trial  pvri0318. 
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Fig.  15.  Range-corrected  s-polarization  lidar  return  within  FOV 
6i+i  as  a  function  of  FOV  0i+1  for  a  distance  of  101  m  (penetration 
depth  of  6  m)  and  an  optical  depth  of  0.4.  the  lidar  return  from 
a  solid  target  located  at  a  distance  of  200  m  is  also  plotted. 
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Fig.  13.  Range-corrected  lidar  return  p-  and  secondary 
s-polarization  lidar  return  as  functions  of  distance  for  0.75-  and 
12-mrad  total  FOVs  for  the  fog-oil  trial  foh0301. 
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Fig.  16.  Range-corrected  s-polarization  lidar  return  within  0£+1 
as  a  function  of  FOV  0i+1  for  distances  of  101  and  105  m  (penetra¬ 
tion  depths  of  6  and  10  m)  and  optical  depths  of  0.2,  0.4,  and  0.8. 


FOV  (mrad) 

Fig.  14.  Range-corrected  s-polarization  lidar  return  within  FOVs 
0;  and  0,+1  as  a  function  of  FOV  0j+1  for  a  distance  of  101  m 
(penetration  depth  of  6  m)  and  an  optical  depth  (O.D.)  of  0.4. 


the  unsealed  lidar  signal  obtained  from  a  solid  target; 
in  this  case  there  is  no  multiple  scattering,  and  the 
displayed  signal  is  the  laser  footprint  with  no  artifact. 

We  use  the  s-polarization  measurements  to  gener¬ 
ate  the  S  vector  of  Eq.  (5),  and  the  matrix  elements 
are  calculated  with  Eq.  (6).  The  particle-size  volume 
density  distributions  q3  are  then  obtained  through 
the  matrix  inversion  of  Eq.  (5)  by  a  conventional 
second-difference  constrained  linear  inversion  tech¬ 
nique.  Figure  17  shows  the  retrieved  q3  for  the  two 
water-droplet  trials  as  well  as  a  typical  Malvern  par¬ 
ticle  sizer  result.  The  results  from  trial  number 
pvh0305  at  optical  depths  of  0.2  and  0.4  match  the 
Malvern  result  fairly  well.  However,  the  results 
from  trial  pvri0318  show  droplets  definitely  smaller 
than  the  Malvern  result.  Figure  18  compares  the 
retrieved  q3  obtained  for  the  fog-oil-droplet  measure¬ 
ments  with  the  result  obtained  from  the  literature.19 
Also  plotted  in  the  same  figure  is  the  retrieved  q3 
from  synthetic  data  generated  with  the  original  PSD 
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Fig.  17.  Comparison,  of  particle-size  density  distributions  ob¬ 
tained  in  two  trials  by  application  of  matrix  inversion  to  secondary- 
polarization  MFOV  measurements  from  water-droplet  clouds  with 
measurements  made  with  the  Malvern  particle  sizer.  O.D.,  opti¬ 
cal  depth. 
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Diameter  ({jm) 

Fig.  18.  Comparison  of  the  particle-size  density  distribution  ob¬ 
tained  by  application  of  matrix  inversion  to  the  secondary- 
polarization  MFOV  measurements  from  fog-oil  cloud  with  a 
published19  distribution. 

from  the  literature.  As  we  pointed  out  in  Section  3, 
the  inversion  technique  does  not  retrieve  the  initial 
distribution  well  for  small  particles;  the  density  dis¬ 
tribution  is  broadened  and  shifted  toward  larger  par¬ 
ticles.  Therefore  it  is  not  surprising  that  the  q3 
retrieved  from  experimental  data  shows  larger  par¬ 
ticles  than  the  values  reported  in  the  literature.  But 


it  is  evident  that  the  measured  distributions  are  sig¬ 
nificantly  wider  than  the  assumed  distribution,  indi¬ 
cating  that  the  actual  oil  droplets  for  the  smoke 
generator  are  larger  than  those  reported  in  the  liter¬ 
ature.  To  quantify  the  agreement- disagreement 
between  the  reported  density  distributions  of  Figs.  17 
and  18  we  list  in  Table  3  the  mean  diameters  calcu¬ 
lated  from  the  retrieved  distributions.  We  trans¬ 
formed  volume  density  distribution  q3  into  number 
density  distribution  q0  to  calculate  the  mean  particle 
diameter  and  the  surface-volume  mean  diameter  (the 
Malvern  and  literature  values  are  given  in  parenthe¬ 
ses).  For  measurements  pvh0305,  the  agreement  is 
fairly  good  between  the  s-polarization  lidar  results 
performed  at  6-  and  10-m  penetration  depths  (optical 
depths  of  0.2  and  0.4)  and  those  obtained  with  the 
Malvern  particle  sizer.  However,  as  we  pointed  out 
above,  measurements  pvri0318  at  the  entrance  of  the 
aerosol  chamber,  i.e.,  at  a  penetration  of  6  m,  lead  to 
smaller  droplets  than  those  obtained  at  the  middle  of 
the  chamber  (lidar  at  105  m  and  the  Malvern  instru¬ 
ment).  This  difference  is  probably  caused  by  the 
evaporation  of  the  water  droplets  in  contact  with  the 
dryer  outdoor  air;  the  relative  outdoor  humidity  went 
down  to  60%  for  this  particular  trial.  The  results 
obtained  for  the  fog  oil  are  in  agreement  within  a 
factor  of  2-3  with  those  obtained  from  the  literature. 
This  is  a  good  result,  as  we  know  that  the  literature 
values  do  not  consider  a  specific  fog-oil  generator  and 
that  the  inversion  technique  is  less  accurate  for  par¬ 
ticles  smaller  than  the  probing  wavelength.  Small 
particles  are  characterized  by  large  scattering  angles, 
and  the  inaccuracy  of  the  inversion  is  possibly  caused 
by  undersampling  of  large  scattering  angles.  For 
example,  less  than  10%  of  the  cloud’s  probed  length 
contributes  to  measurements  at  scattering  angles 
greater  than  45°  for  our  measurements  made  at  a  6-m 
penetration  depth. 


5.  Limitation  of  Present  Model 

In  Section  1  we  pointed  out  that  our  model  is  re¬ 
stricted  to  small  optical  depths,  so  first-order  forward 
scattering  is  predominant.  To  quantify  the  small 
optical  depth  restriction,  we  estimate  the  relative 
contributions  of  the  different  scattering  orders  as 
functions  of  optical  depth.  We  proceed  as  in  Ref.  1 
and  use  an  extension  of  the  Beer— Lambert  law. 
Given  I0,  the  incident  intensity  on  a  homogeneous 
cloud,  the  unscattered  intensity  loss  and  the  amount 


Table  3.  Comparison  of  the  Measured  Mean  Values® 


Trial  Number 

Penetration 
Depth  (m) 

Optical  Depth 

{x)qs  (|xm) 

(x)q0  (p,m) 

(xs)q0/{x2)q0  (pm) 

(Pm) 

pvri0318 

6 

0.4 

14.4  (28.8) 

6.9  (8.2) 

11.7  (18.9) 

14.9  (17.3) 

pvh0305 

6 

0.2 

22.5  (28.8) 

8.5  (8.2) 

16.7  (18.9) 

15  (17.3) 

pvh0305 

10 

0.4 

23.8  (28.8) 

11.6(8.2) 

19.6  (18.9) 

22.6  (17.3) 

fori0308 

6 

0.4 

2.5  (1) 

0.3  (0.3) 

1.4  (0.67) 

1.2  (0.8) 

foh0301 

6 

0.4 

2.4  (1) 

0.3  (0.3) 

1.4  (0.67) 

1.2  (0.8) 

“The  values  in  parentheses  were  obtained  with  the  Malvern  Particle  Sizer  or  from  the  literature. 
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of  scattered  intensity  as  functions  of  penetration 
depth  obey  the  Beer-Lambert  law,  i.e., 


It  =  exp(-ctz),  (10) 

Idi  =  fo[l  ~  exp(-ocz)],  (11) 

where  Ir  is  the  transmitted  intensity  and  ID1  is  the 
scattered  intensity.  Part  of  the  scattered  intensity 
is  rescattered  and  the  process  is  repeated  for  the 
higher  scattering  orders .  In  the  small-angle  scatter¬ 
ing  approximation  (predominance  of  forward  scatter¬ 
ing),  the  transmission  loss  of  the  doubly  scattered 
beam  can  be  approximated  by 

AItdi  =  -/yj5i[l  -  exp(-az)]aAz,  (12) 

and  the  solution  is 

Itdi  =  Iq  exp(-az){exp[l  -  exp(-ctz)]},  (13) 
Id2  =  /o[1  -  exp(-oz){exp[l  -  exp(-az)]},  (14) 

where  Itdi  is  the  fraction  of  the  scattered  beam  (ID1) 
that  is  transmitted  and  ID2  is  the  fraction  that  is 
rediffused.  Proceeding  in  the  same  way  to  third  or¬ 
der,  we  find  that 

Itd2  =  h  exp(-«2  -  l)exp{exp[l  -  exp(— ctz)]},  (15) 

Id3  —  /o(l  ~  exp (—ae  -  l)exp{exp[l  -  exp(-ctz)]}), 

(16) 

where  ITD2  is  the  fraction  of  the  scattered  intensity 
ID2  that  is  transmitted  and  Im  is  the  fraction  that  is 
rediffused.  Finally,  the  transmission  loss  to  fourth 
order  is  written  as 


ItD3  =  /o  exP 

Id4  =  /o(l  ~  IrTDa)>  (18) 

where  ITD3  is  the  fraction  of  the  scattered  intensity 
Id 3  that  is  transmitted  and  ID4  is  the  fraction  that  is 
rediffused. 

Figure  19  shows  the  relative  contributions  of  the 


-  IDsa(z) 


dz, 


(17) 
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Fig.  19.  Approximate  relative  lidar  contributions  of  the  three 
scattering  orders  as  functions  of  apparent  optical  depth. 


various  scattering  orders  (ID4-i  -  ID  l)  normalized 
with  scattered  intensity  ID1  as  functions  of  the  optical 
depth  divided  by  2.  These  expressions  model  the 
relative  contributions  by  higher-order  forward  scat¬ 
tering  to  lidar  returns  in  the  approximation  that 
backscatter  sensitivity  to  scattering  order  is  slight. 
The  division  by  2  is  necessary  because  in  lidar  geom¬ 
etry  the  two-way  transmission  through  the  scattering 
medium  has  to  be  considered.  The  higher  scattering 
orders  successively  become  predominant  as  optical 
depth  increases.  Clearly,  the  contributions  of 
second-order  forward  scattering  cannot  be  neglected 
for  optical  depths  greater  than  0.3-0.4.  The  effects 
of  higher  forward-scattering  orders  are  to  spread  the 
scattered  photons  at  larger  angles  and  also  to  redirect 
energy  back  into  the  inner  FOY^,  making  the  rela¬ 
tionship  given  by  Eq.  (1)  between  the  scattering  angle 
and  the  FOV  invalid.  How  is  this  extra  energy  in¬ 
terpreted  by  the  present  model?  Can  the  higher 
scattering  orders  be  properly  modeled?  These  are 
questions  outside  the  scope  of  this  paper,  but  they 
must  be  answered  if  the  PSD  retrieval  from  lidar 
secondary  polarization  is  to  be  expanded  to  optical 
depths  greater  than  0.4. 

The  particle-size  resolution  is  range  and  size  de¬ 
pendent.  The  resolution  of  large  particles  requires 
small  FOV7 s,  and  those  FOVs  have  to  be  even  smaller 
when  the  cloud  is  situated  at  a  greater  distance  from 
the  lidar.  Small  FOV5 s  lead  to  the  following  require¬ 
ments  or  restrictions: 

The  FOV  ring-iris  mask  disk  must  be  positioned  in 
the  image  plane  of  the  principal  optics.  To  cover 
cloud  near  and  far  distances  and  maintain  FOV  pre¬ 
cision  it  is  necessary  to  adjust  the  image  plane  posi¬ 
tion  with  a  translation  stage. 

Ideally  the  laser  spot  size  on  the  cloud  should  sub¬ 
tend  an  angle  smaller  than  the  smallest  FOV.  This 
means  that  outgoing  laser  beam  size  and  divergence 
smaller  than  25  mm  and  0.1  mrad,  respectively,  are 
required. 

The  transmitter  and  receiver  optics  has  to  be  col- 
linear,  and  the  alignment  must  be  maintained  within 
the  smallest  FOV. 

6.  Conclusion 

The  proposed  technique  based  on  the  use  of  the  sec¬ 
ondary  polarization  offers  good  potential  for  the  de¬ 
termination  of  the  PSD  of  clouds.  (The  advantage  is 
that  the  useful  multiple-scattering  signal  is  not  bur¬ 
ied  under  strong  single-scattering  contributions. 
The  proposed  method  is  applicable  at  small  penetra¬ 
tion  depths,  but  this  is  the  very  region  where  it  is 
useful,  i.e.,  where  it  is  difficult  to  separate  multiple 
from  single-scattering  contributions  in  conventional 
total  return  measurements.  The  present  method  is 
applicable  only  at  small  penetration  depths  and  does 
not  require  absolute  knowledge  of  extinction  coeffi¬ 
cient  a(z),  which  is  set  proportional  to  the  measured 
total  lidar  signal.  The  model  is  relatively  simple  in 
application,  and  the  measurements  are  easily  done 
with  FOV  segmentation.  The  measurements  re¬ 
quire  low  laser  beam  divergence,  small  laser  beam 
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size,  and  imaging-quality  collection  optics.  The 
technique  can  be  applied  for  particle  sizing  of  large 
spray-nozzle  systems,  stack  emissions,  and  thin  nat¬ 
ural  clouds.  Finally,  because  the  technique  uses  a 
unique  characteristic  of  spherical  droplets,  i.e.,  zero 
depolarization  at  180°,  it  cannot  be  readily  applied  to 
nonspherical  particles  or  to  large  droplets  that  are 
frequently  found  to  be  oblate. 
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